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NOMENCLATURE 
Surface area of each element 
Expansion coefficient 
gravity 
Thermal conductivity 
Area divided by the perimeter 
Rate of energy output from each element 
Net rate of energy input from each element to the fluid 
Lead resistance 
Precision resistor resistance 
Ambient temperature 
temperature of the element 
Voltage measured over the leads and heater [Volts] 


Voltage over the precision resistor, leads and heater [Volts] 


Kinematic viscosity 
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I. INTRODUCTION 


As electronic components become increasingly smaller, the volumetric 
power generation rates escalate and the problem of their thermal management 
becomes even more challenging. The greatest constraint to further advances 
in decreasing the feature size and increasing functional capability of electronic 
components could be the unavailability of appropriate heat removal techniques. 
The performance and reliability of microelectronic components is a function of 
junction temperature. The average failure rates of microelectronic components 
increase exponentially with junction temperature increase; Bar-Cohen and 
Kraus [Ref. 1]. 

Direct liquid cooling has emerged as one of the most promising thermal 
management techniques for overcoming this barrier; Bar-Cohen [Ref. 2]. The 
advantages of liquid cooling have been demonstrated for over 40 years in 
macroelectronic components, yet its use for cooling microelectronic components 
is relatively new. 

There are three types of liquid cooling techniques that can be employed; 
boiling, forced convection, and natural convention. Initiation of boiling with 
dielectric fluids is complicated by the highly wetting nature of these liquids. 
Forced convection can produce higher heat dissipation rates, however, it may 


present its own set of problems with reliability and higher costs. Even though 


natural convection cooling typically allows lower heat dissipation rates than 
the first two methods, it has the advantages of simplicity of design, low 
operation cost due to no outside power requirements and minimum 
maintenance, absence of noise and high reliability. Also, in the event of a 
mechanical failure, natural convection may be the only method of cooling 
available [Refs. 3, 4, 5]. 

Natural convection originates when a body force acts on a fluid in which 
there are density gradients. The net effect is a buoyancy force, which induces 
free convection currents. In the most common case, the density gradient is due 


to a temperature gradient, and the body force is due to the gravitational field. 


A. SUMMARY OF LITERATURE 

Baker [Ref. 6] studied the heat transfer characteristics of flush mounted 
heat sources on vertical surface in several liquids. Park and Bergles [Ref. 7] 
investigated both flush and protruding arrangements in water and R-113. Bar- 
Cohen and Schweitzer [Ref. 8] presented correlations for natural convection in 
a vertical channel dissipating uniform heat flux from the sides. Kelleher ef. al. 
[Ref. 9] conducted flow visualization and heat transfer measurements for 
natural convection from a long heated protrusion on a vertical insulated wall 
of a rectangular chamber filled with water. Lee et. al. [Ref 10] supported 


Kelleher’s findings with a computational investigation. 


Joshi et. al. [Ref. 11] presented flow visualizations and component surface 
temperature measurements for natural convection cooling of a three by three 
array of discrete heated protrusions on the vertical wall of a rectangular 
enclosure filled with a dielectric liquid. Joshi et. al. presented experimental 
studies of the heat transfer and flow characteristics of a column of protruding 
heat sources, [Ref. 12] on a vertical surface and within a vertical channel, [Ref. 
13]. Sathe and Joshi [Ref. 4] reported results of a numerical investigation of 
natural convection flow and heat transfer arising from a protruding heat 
source on a vertical plate within an enclosure. Joshi and Paje [Ref. 14] 
reported experimental results of natural convection heat transfer from a 
commercially available semiconductor device package. 

Joshi and Knight [Ref. 15] studied the natural convection from a single 
column of eight in-line, rectangular heat sources flush mounted on one wall of 
a vertical channel immersed in water. Gaiser [Ref. 16] examined natural 
convection adjacent to a three column array of 15 per column uniformly heated 
components flush mounted on a vertical flat plate in water. Haukenes [Ref. 5] 
reported results of a numerical investigation of two dimensional natural 
convection flow and heat transfer from a substrate-mounted flush heat source 
immersed in a liquid filled square enclosure. Akdeniz [Ref. 17] investigated the 
natural convection heat transfer from an array of heaters flush mounted on a 


vertical test surface in response to both step and periodic input power. 


B. PRESENT STUDY 

The investigation reported here is a continuation of the studies of Gaiser 
[Ref. 16], Haukenes [Ref. 5], and Akdeniz [Ref. 17]. As stated above, Akdeniz 
[Ref. 17] presented the changes in the natural convection flow and heat 
transfer characteristics of a three column array of 45 flush mounted heaters 
on a vertical test surface. In his study, power to the bottom heater of the 
center column oscillated with time while power to all the other heaters was 
kept constant. The present experimental study was conducted with the same 
test surface as in [Ref. 16, 5, and 17], however, the entire center column was 
powered in a time periodic fashion while the side columns of heaters remained 
unpowered. 

Specifically, the natural convection heat transfer response of a column of 
heaters flush mounted on a vertical test surface in water to periodic input 
power has been investigated. Two types of periodic input power variations were 
examined: a triangular wave and an approximate square wave. The resulting 
heater temperatures over several cycles were measured for mean values of 0.5, 
1.0, 2.0, and 3.0 watts and varying amplitudes. The frequency of the input 
power patterns was also varied, from 0.025 to 0.1 Hz. The measured heater 
temperatures were compared with the responses for steady input power equal 


to the mean of the periodic input. 


C. OBJECTIVES 
The particular objectives of this investigation were: 
¢ To examine the effects of power pulsation on the heat transfer 
characteristics for the selected conditions. 


¢ To determine if any enhancement of heat transfer occurs due to the input 
power variations. 


Il. EXPERIMENT 


A. EXPERIMENTAL APPARATUS 

The experimental apparatus is described in detail by Gaiser [Ref. 16]. 
Additions and modifications to the apparatus were discussed by Haukenes 
(Ref. 5] and Akdeniz [Ref. 17]. The following is a brief summary of the above 
with emphasis on the additional changes to reflect the current status of the 
apparatus. 

The apparatus itself is sub-divided into four separate assemblies; test 
surface assembly, power supply assembly, data acquisition/reduction unit, and 
the deaerating assembly, as seen in Figure 1. The first three of these systems 
are interrelated through three distinct variables; the input power to the heater 
strips, the heater strip surface temperatures, and the ambient temperature of 
the bath. The last of these systems, the deaerating system is used separately 


to reduce the amount of entrained air and compressible gasses. 


1. Test Surface Assembly 
This consisted of a vertical test surface with 45 flush mounted 
heater strips placed on three columns of 15 each. The test surface is immersed 
in a one cubic meter plate glass tank filled with deionized water. The water is 
deionized by a series of filters attached externally to the tank. The test surface 


was constructed using a 31.75 cm x 31.75 cm x 1.11 cm plexiglass plate with 
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Figure 1. Experimental Apparatus. 


3.58 cm x 0.79 cm slots milled out to a depth of 0.19 cm at each location where 
a heater strip was to be embedded. Figure 2 shows the front view of the test 
surface showing the location of each heater strip in the central column, 
powered in this investigation. Within each slot, a 0.79 cm diameter hole was 
drilled through the board in order to allow the heater strip power leads and 
thermocouple wires to pass through. Side and top views of a slot are shown in 
Figure 3. [Ref. 16] 

Before mounting the heater strips to the board, two power leads and 
a 0.0076 cm copper-constantan thermocouple were attached to the back of each 
strip. The power leads were simply soldered onto the heater strip. A layer of 
Omega bond 101 high thermal conductivity adhesive was coated over the back 
of the heater strip to electrically isolate each thermocouple from the heater 
strip, while allowing for accurate temperature measurements. Once the strip 
was coated, an additional small drop of the adhesive was used to attach the 


thermocouple onto the center of the heater strip. [Ref.16] 


2. Power Supply Assembly 
The power supply assembly provided power to the individual heater 
strips. For this investigation only the central column was used as seen in 
Figure 2. An HP-85 computer was used in conjunction with an HP-59501 
power supply programmer and a 0-40 Volt, 0-1.56 Amp HP-6298A D.C. power 
supply to control the power output of the center column of heaters. The HP- 


6298A power supply was connected to two of the six distribution panels. These 
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Figure 2. Front View of Test Surfaces with Active Heaters. All 
Units in cm. 








Figure 3. Front and Side View of the Test Surface Heater Slots. 
All Units in cm. 
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two distribution panels are the ones to which the center column of heaters are 
connected. Each distribution panel acts as an individual power supply for a 
specified number of heaters. Each panel supplies power to either seven or eight 
heaters. If only seven heaters were connected to a particular panel, a ten ohm 
resistor was used to equilibrate the total power requirement of the panel as 
compared with an eight heater panel. 

Figure 4 shows the power distribution panel. Each heater on a panel 
was connected in series with a 2.01 ohm 1% precision resistor. The heater and 
resistor pair was then connected in parallel to two copper wire bus bars, 
previously referred to as distribution panels. Additional 1% and 5% precision 
resistors were connected in parallel with the original 1% precision resistors in 
order to approximate a constant input power for all heaters. 

. The data acquisition system measures the total voltage drop across 
the bus, V,, and the voltage drop for each heater including the leads, V,. Using 
these two measured quantities the power output for each heater was 


determined using the following formula: 


V(R._-R,)- V 
Power = Sate Taha Z v i (V> - 


) 
ft 
Toy 


Where R, is the precision resistor resistance and R, is the heater resistance 


plus lead resistance. [Ref. 5] 
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3. Acquisition/Reduction Unit 

The data acquisition/reduction unit was made up of an HP 3852A 
Data Acquisition unit, an HP 9153 Computer system, and an IBM clone 
personal computer. This system initiates, records and processes all data from 
the test surface and power assembly. When tasked by the HP 9153 Computer 
system, the HP 3852A Data Acquisition unit reads the thermocouple voltages, 
and the appropriate voltages for the determination of heater powers. 

This data was then transmitted into the HP 9153 Computer system 
which processed it into temperatures and power measurements. The 
temperature readings were further manipulated by taking the difference 
between the ambient temperature and the temperature of the heater. The 
power readings were processed as per the discussion above. 

This data, once processed, was transferred to an IBM clone personal 
computer. This was accomplished by printing the data directly to the RS-232 
port of the HP 9153 Computer system, which was connected to the RS-232 
IBM computer via a cable. The statement used for the HP 9153 computer to 
transmit through the RS-232 port was "printer is 9". The HP 9153 computer 
automatically converted the data file into ASCII file format prior to printing 
it to the RS-232 port. The data was received into the IBM computer by using 
a communication program. Once received by the IBM computer, the data was 
imported into a spread sheet program which was used for graphics. The data 


files were imported as a comma & "" delimited file. The data file had to be 


13 


formatted to be imported this way. This was accomplished using an editor 


within the IBM clone computer. 


4. Deaerating Assembly 

The deaerating assembly was composed of a vacuum pump, 
discharge pump, a heating element, and the deaerating tank, see Figure 5. The 
deaerating tank can be used in conjunction with the filtration system. The 
filtration system, consisting of four cartridges, is used to ensure the purity by 
maintaining the resistivity of the tank water at 0.7 megaohm-cm or above. 
This is achieved through periodic recirculation of the tank water. 

The bath water from the tank was vacuum dragged into the 
deaerating tank through a nozzle. The vacuum was manually regulated to be 
between 68 and 85 KPa. A heating element within the deaerating tank heats 
the water to 10 - 15 K above the bulk temperature in the test tank. The 
temperature within the deaerating tank is maintained by regulating the flow 
rate through the tank by the discharge pump. Flow rates between 900 
milliliters per minute to 1500 milliliters per minute were the norm. 

Float switches were used to maintain the water level within a range 
around the nominal level. The capacity of the deaerating tank was 60.6 liters 
and the tank was nominally maintained to be half full. This was done to 
provide complete coverage of the panel heating element. The tank itself was 
elevated 1.83 meters above the floor to provide sufficient head for the discharge 
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B. EXPERIMENTAL PROCEDURE 

In preparing for a series of experiments, the mechanical stirrer would be 
run for five to ten minutes to dissipate prior stratification within the test tank. 
Crushed ice was placed into the thermocouple reference bath. Typically the 
data acquisition unit and the HP computers were left on. Only the IBM 
computer and the power supply were secured. Quiescence was achieved within 
the test tank prior to powering any of the heaters. This would take 
approximately 30 minutes. 

Once the bath was ready, the power supply was energized and the 
program to generate appropriate periodic power input patterns through the 
heaters was loaded and run (see Appendix D). Another 10 to 12 minutes were 
needed to ensure that steady state conditions were reached prior to taking any 
temperature measurements. Once steady conditions were achieved, the data 
acquisition program was started and periodic temperature measurements over 
five minutes were taken. 

The effects of two different type of periodic input power patterns were 
examined: a triangular pattern and an approximate square wave pattern. 
Measurements at three locations along the test surface were obtained for each 
run. These measurements were obtained for four mean input power levels, 


each for three amplitudes and three frequencies of pulsation. The 
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corresponding steady state measurements for the selected mean power levels 


were also made. 


iW 


III. EXPERIMENTAL RESULTS 


A. DESCRIPTION OF DATA 

As previously mentioned, the test surface consisted of three columns of 
15 heaters, numbered as in Figure 6. Only the center column of heaters was 
used for this study, heaters 16 through 30. Three different configurations of 
the center column of heaters were used, see Figure 7. In configuration A, all 
15 heaters were powered uniformly. In configuration B, only the top 13 heaters 
were powered uniformly, thus leaving heater 29 and 30 de-energized. 
Configuration C is a further modification of heater configuration B, de- 
energizing the top heater, heater 16. This left 12 heaters heated. The reason 
for the modifications was to the failure of some heaters due to prolonged 
periodic input power through them. 

Contained within Appendix A is a complete listing of all runs performed 
for this study. All graphical data is presented within Appendix B. The captions 
for each figure contained within Appendix B refer to a heater configuration, 
using a three digit alphanumeric format. The first digit is a letter, either A, B, 
or C, which references the number of active heaters as stated above, see Figure 
7. The second and third digits are numbers which denote the heater number, 


as seen in Figure 6. For example, heater configuration B21 means that the top 
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Figure 7. Heater Configurations A, B, and C. 
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13 heaters where heated and the heater of interest is heater number 21, or the 
middle heater. 

For this investigation, measurements at three heater locations were made 
for each run near the top, middle and bottom of the column, in order to 
determine the downstream effects on the heat characteristics. Table I, shows 
which heaters were used near the top, middle, and, bottom locations for each 


of the heater configurations. 


TABLE I. DESIGNATION OF TOP, MIDDLE AND BOTTOM 
HEATERS FOR DIFFERENT HEATER 
oe 







HEATER HEATER CONFIGURATION 


MIDDLE 


a 


Each of the figures found in Appendix B are numbered with a two letter 
code. The first letter denotes a specific mean power input and amplitude. The 
second letter designates a specific heater, top, middle, or bottom, and frequency 
at which the input power pattern was applied. This can be seen in Table II and 
Table III. It takes a grouping of three figures to represent a single run. A run 
represents a grouping of figures where only the heater of interest changes for 
a given set of power input parameters. The second letter designation for the 


figures, A, B, and C, constitutes one run. E, D, and F constitutes the next run, 


zal 


and G, H, and I| constitutes the last run. These groupings are divided up by the 
frequency at which the power pattern is applied. 
TABLE II. FIRST LETTER ine OF FIGURE CODE 


First Letter|. Mean, misAmphitude |, PowemPatieaal 
Designation Value of 
Power 
Input 


A 
PB | 2 set} Tritngulanst 
ee 
oe 
ee 

PG | sd | Triangular 
Pom | st | Soianguter 
a i PE bass 

Square 

a 
[os [Square 


Square 
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To obtain all the data for a single run, each experiment had to be 


repeated a total of six different times. This allowed gathering the necessary 
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TABLE III. SECOND LETTER DESIGNATION FOR FIGURE CODE 


Second Letter Heater Frequency 
Designation 


; 0.100 
B 0.100 
c 0.100 

0.050 


0.050 
0.050 
Top 0.025 

Middle —«| 0.025 





input power and temperature data on three different heaters with pulsating 
and steady input powers. The pulsating and steady programs for the HP 85 
computer are contained within Appendix D. The data are presented one heater 
per graph. This was done to enhance readability. 

For all discussion of the data two definitions are needed. The first being 
heat transfer enhancement; where the mean cyclic temperature response is 
less than the steady temperature response. The second being significant heat 
transfer enhancement; where the maximum cyclic temperature response is less 


than the steady temperature response. 
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B. DISCUSSION OF TRIANGULAR WAVE DATA 

The baseline run for the triangular wave data was arbitrarily taken to be; 
mean value power input of 2.0 watts, amplitude of 1.0 watts, and a frequency 
of 0.050 Hz, see Figures BD-BF. 

The amplitude of the cyclic temperature response as a percentage of the 
steady temperature response ranges from 6.7% at the top heater to 16.8% at 
the bottom heater, where the amplitude of input power is 50% of the mean 
power input. There are some damping effects evident from the temperature 
response. The damping effects are more pronounced at the top heater than the 
bottom or middle heaters. There seems to be little change in the damping 
effects between the middle heater and the bottom heater. Minimizing the 
amplitude of the temperature response is desirable for application purposes. 
The frequency of the cyclic temperature response matches that of the power 
input, however, as expected there was a phase lag. Also, for the top heater the 
cyclic temperature response, is slightly lower than the steady temperature 
response implying some heat transfer enhancement. The pulsation of the power 
input seems to have a more dramatic effect on the top heater than the bottom 
heater, where there seems to be little to no heat transfer enhancement. This 


trend will become more dramatic later on. 


1. Effects of the Amplitude of Power Pulsation 
Figures CD-CF present results for mean input of 2.0 watts, 


amplitude of 2.0 watts and a frequency of 0.025 Hz. These figures are an 
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increase in the amplitude of 1.0 watts compared to the baseline conditions. For 
these conditions the heat transfer enhancement is more pronounced than the 
base line case. For the top heater the mean cyclic temperature response, is 6 
K lower than the steady temperature response. The top heater sees the 
greatest enhancement and the bottom heater the least. For the top and middle 
heaters, the maximum cyclic temperature response is clearly lower than the 
steady temperature response. Note, the maximum cyclic temperature response 
for the top and middle heaters are very close. For the bottom heater, the mean 
of the cyclic temperature response is significantly lower than the steady 
temperature response. Still, the maximum cyclic temperature response is 
higher than the steady temperature response. 

Figures AD-AF present results for mean input of 2.0 watts, 
amplitude of 0.5 waits and a frequency of 0.050 Hz. These figures are a 
decrease in the amplitude of 0.5 watts compared to the baseline case. The heat 
transfer enhancement effects due to this pulsation are negligible. The 
amplitude of the temperature response is less than in the baseline case. 

It can be readily seen that as the amplitude to mean power ratio 
approaches one, significant heat transfer enhancement takes place. The 
importance of this ratio will be demonstrated again later on with different 


input powers. 
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2. Effects of the Frequency of Pulsation 

Figures BA-BC present results for mean input of 2.0 watts, 
amplitude of 1.0 watts and a frequency of 0.100 Hz. These figures are an 
increase in the frequency by a factor of two from the baseline case. The first 
effect noticed is the damping of the temperature response. The amplitude 
response of the top heater drops from 6.7% in the base line case to 3.7%, and 
bottom heater drops from 16.8% to 6.4%. Some increase in enhancement is 
noted by increasing the frequency. Figure BA shows a slightly lower mean 
cyclic temperature response than the steady temperature response. 

Figures BG-BI present a decrease in the frequency of the pulsation 
of the input power from 0.050 to 0.025 Hz. The effects of damping were 
diminished. The temperature amplitude at the top heater jumped from 6.7% 
to 10% and for the bottom heater from 16.8% to 20.5%. With the decreased 
frequency the thermal inertia of the surface and thermocouple becomes less 
important. No significant heat transfer enhancement was noted at the bottom 
heater location. A slight enhancement at the top heater was noted, but less 
than in the base line case. 

In summary, as the frequency of the power pulsation is increased the 
amplitude of the temperature response decreases and there is a slight heat 
transfer enhancement. As stated previously, smaller amplitudes of the 


temperature response are desirable for applications purposes. 
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3. Effects of the Mean Input Power 

Figures HD-HF present results for a mean input of 3.0 watts, 
amplitude 1.0 watts and a frequency of 0.050 Hz. These figures are an increase 
in the mean input power of 1.0 watts from the baseline case. It is readily seen 
that the amplitude of the temperature response remains unchanged from the 
baseline case for the top and bottom heater locations. For the middle heater 
location the amplitude of the temperature response is slightly larger. This may 
be due to the fact that heater configuration B for this run was used instead of 
heater configuration A for the baseline run resulting in minor differences in 
the response times of the embedded surface thermocouple. The frequency 
response still matches that of the input power for all three heaters. The drop 
in temperature 1s more pronounced for the top heater than in the base line 
case. 

Figures FD-FF present results for a mean input of 1.0 watts, 
amplitude of 1.0 watts and a frequency of 0.050 Hz. Heat transfer 
enhancement is seen for all three heaters. The greatest enhancement occurred 
for the top heater location. Note heater configuration B was used for this run. 
The frequency response matches that of the input power, and temperature 
amplitudes are approximately the same as that of the baseline case. The 
middle heater location has a slightly larger temperature amplitude than in the 
baseline case, probably due to the different heater configuration. For the top 


heater location there is approximately a 4 K difference between the mean of 
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the cyclic temperature response and the steady temperature response, see 
Figure FD. 

The effects of the mean power are two fold. First as a ratio of the 
amplitude to mean power level. Also, the mean power level is important as an 
independent variable. The larger the mean power input the more vigorous the 
base flow is which is more prone to amplification of disturbances. These 
influences can be seen in Figures FD-FF and Figures HD-HF. In Figures FD- 
FF the amplitude to mean power ratio is equal to one. Significant heat transfer 
enhancement is occurring. In Figures HD-HF, the amplitude to mean power 
input ratio is equal to 0.33, still greater heat transfer enhancement can be 
seen than in the baseline case where the amplitude to mean power ratio is 
equal to 0.50. This is accounted for by the higher mean input power than the 


baseline case. 


4. Combined of Effects 
Itis also important to examine more than one variation to the input 
power pulsation, (e.g. increase the frequency and the amplitude of the input 
power keeping the mean input the same). 
a. Varying Amplitude and Frequency 
Figures CA-CC present results for mean input of 2.0 watts, 
amplitude of 2.0 watts and a frequency of 0.100 Hz. These figures are an 


increase in the amplitude of 1.0 watts and an increase of the input frequency 
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by a factor of two from the baseline case. There was significant heat transfer 
enhancement for all three heater locations. Significant heat transfer 
enhancement was previously defined as the maximum cyclic temperature 
response being less than the steady temperature response. The amplitude of 
the temperature response is greater than the baseline case, however, less than 
if the amplitude of the input power alone was increased, see Figures CD-CF. 

Figures CG-CI present results for mean input of 2.0 watts, 
amplitude of 2.0 watts and a frequency of 0.025 Hz. These figures are an 
increase in the amplitude of 1.0 watts and a decrease of the input frequency 
by a factor of two from the baseline case. Heat transfer enhancement was 
observed for all three heater locations. For the middle and bottom heater 
locations the mean cyclic temperature response was well below the steady 
temperature response. As expected, the amplitude of the cyclic temperature 
response was significantly larger than the baseline case and greater than if the 
frequency of the pulsation alone was decreased, see Figures BG-BI. 

Figures AA-AC present results for mean input of 2.0 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These figures are a 
decrease in the amplitude of 0.5 watts and an increase of the input frequency 
by a factor of two from the baseline case. Heat transfer enhancement was 
found for the top and bottom heater locations. The maximum cyclic 
temperature response for the top heater location coincided with the steady 


temperature response. For the middle heater the mean cyclic temperature 


Vés) 


response was equal to the steady temperature response indicating no heat 
transfer enhancement. For the bottom heater location the mean cyclic 
temperature response was slightly lower than the steady temperature 
response. The amplitude of the temperature response was significantly smaller 
than the baseline case for all three heater locations tested. 

Figures AG-AI present results for mean input of 2.0 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These figures are a 
decrease in the amplitude of 0.5 watts and a decrease of the input frequency 
by a factor of two from the baseline case. No heat transfer enhancement was 
found for any of the three heater locations. The amplitude of the temperature 
response was slightly less than the baseline case. 

The variation of amplitude and frequency act independently of 
each other. In fact, the benefits of higher amplitudes and higher frequencies 


to the heat transfer appear to be additive. 


6. Varying Amplitude and Mean Input Power 
Figures ID-IF present results for mean input of 3.0 watts, 
amplitude of 2.0 watts and a frequency of 0.050 Hz. These figures are an 
increase in the amplitude of 1.0 watts and an increase of the mean power of 
1.0 watts from the baseline case. Heat transfer enhancement is seen for all 
three heater locations and significant heat transfer enhancement occurs for the 
top heater location. As expected, the amplitude of the temperature response is 


greater than the baseline case. 
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Figures GD-GF present results for mean input of 3.0 watts, 
amplitude of 0.5 watts and a frequency of 0.050 Hz. These figures are a 
decrease in the amplitude of 0.5 watts and an increase of the mean power of 
1.0 watts from the baseline case. The top and middle heater locations 
demonstrate some heat transfer enhancement, more so for the top heater. The 
bottom heater shows no heat transfer enhancement and even a slight decrease 
in the heat transfer. The amplitude of the temperature response is less than 
in the baseline case. The amplitude of the temperature is greatest for the 
middle heater, of the three heater locations. Note heater configuration B was 
used in this case and typically smaller amplitudes of the temperature response 
for the bottom heater location were found than for heater configuration A. 

Figures ED-EF present results for mean input of 1.0 watts, 
amplitude of 0.5 watts and a frequency of 0.050 Hz. These figures are a 
decrease in the amplitude of 0.5 watts and a decrease of the mean power of 1.0 
watts from the baseline case. The heat transfer enhancement at the top heater 
location is greater than for the baseline case. Here the maximum cyclic 
temperature response is equal to the steady temperature response. While in 
the baseline case it was greater than the steady temperature response. The 
temperature response for the middle and bottom heater locations is 
approximately the same as that of the baseline case. As expected, the 


amplitude of the temperature response was less than for the baseline case. 
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Figures DD-DF present results for mean input of 0.5 watts, 
amplitude of 0.5 watts and a frequency of 0.050 Hz. These figures are a 
decrease in the amplitude of 1.5 watts and a decrease of the mean power of 0.5 
watts from the baseline case. For all three heater locations heat transfer 
enhancement was seen. Significant heat transfer enhancement was seen for 
the top heater location. Once again, this emphasizes the importance of the 
increase in the amplitude to mean power ratio for heat transfer enhancement. 
The amplitude of the temperature, as expected, is less than in the baseline 


Case. 


c. Varying Frequency and Mean Input Power 

Figures HA-HC present results for mean input of 3.0 watts, 
amplitude of 1.0 watts and a frequency of 0.100 Hz. These conditions are an 
increase in the mean power input of 1.0 watts and an increase of the input 
frequency by a factor of two from the baseline case. There was heat transfer 
enhancement for all three heater locations with significant enhancement for 
the top heater location. The amplitude of the temperature response was less 
than that of the baseline case. 

Figures HG-HI present results for mean input of 3.0 watts, 
amplitude of 1.0 watts and a frequency of 0.025 Hz. These conditions are an 
increase in the mean power input of 1.0 watts and a decrease of the input 
frequency by a factor of two from the baseline case. Only the top and middle 


heater locations presented a slight heat transfer enhancement. As expected, 
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the amplitude of the temperature response was greater than that of the 
baseline case. 

Figures FA-FC present results for mean input of 1.0 watts, 
amplitude of 1.0 watts and a frequency of 0.100 Hz. These figures are a 
decrease in the mean power input of 1.0 watts and an increase of the input 
frequency by a factor of two from the baseline case. Significant heat transfer 
enhancement was observed for all three heater locations. The amplitude of the 
temperature response is less than that of the baseline case. This run 
represents dramatic improvements to the heat transfer characteristics using 
lower mean power input and variations to the amplitude. Being able to use 
lower amplitudes and mean power inputs to achieve significant heat transfer 
enhancement, results in lower mean cyclic temperature responses. This type 
of enhancement to the heat transfer characteristics with lower input powers 
would be desirable for applications purposes, when such enhancement could 
be achieved by using auxiliary heaters with pulsed input power. 

Figures FG-FI present results for mean input of 1.0 watts, 
amplitude of 1.0 watts and a frequency of 0.025 Hz. These figures display the 
effect of a decrease in the mean power input of 1.0 watts and a decrease of the 
input frequency by a factor of two from the baseline case. Heat transfer 
enhancement at all three heater locations is evident with significant 
enhancement for the top heater location. Amplitude of the temperature 


response was greater than the baseline case. 
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The frequency of the pulsating power has less of an effect on the 
heat transfer enhancement than the amplitude and mean power input does. 


Its greatest influence is on the amplitude of the temperature response. 


d. Other Variations of Parameters 

Figures DA-DC present results for mean input of 0.5 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These parameters are a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.5 
watts, and an increase in the input frequency by a factor of two from the 
baseline case. There is significant heat transfer enhancement for all three 
heater locations. This again demonstrates the importance of the amplitude to 
mean power input ratio. As expected, the amplitude of the temperature 
response is less than in the baseline case. In addition, for the top heater 
location the mean cyclic temperature response is 20 K below the mean cyclic 
temperature response for the baseline case. The corresponding decrease 18.5 
K and 14 K for the middle and bottom heater locations respectively. This is a 
significantly lower mean temperature response with greater heat transfer 
enhancement than the baseline case. Also, the mean cyclic temperature 
response for the middle and bottom heater locations are the same. 

Figures EA-EC present results for mean input of 1.0 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These values are a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.0 


watts, and an increase in the input frequency by a factor of two from the 
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baseline case. There is heat transfer enhancement for all three heater 
locations. Significant heat transfer enhancement occurs for the top heater 
location only. The amplitude of the temperature response is less than the 
baseline case. 

Figures IA-IC present results for mean input of 3.0 watts, 
amplitude of 2.0 watts and a frequency of 0.100 Hz. These figures are an 
increase in the amplitude of 1.0 watts, an increase of the mean power of 1.0 
watts, and an increase in the input frequency by a factor of two from the 
baseline case. There is heat transfer enhancement for all three heater locations 
with significant heat transfer enhancement for the top heater only. The 
amplitude of the temperature response is approximately the same as the 
baseline case. The reduced amplitude of the temperature response is a function 
of the higher frequency of the input power. 

Figures GA-GC present results for mean input of 3.0 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These figures are a 
decrease in the amplitude of 0.5 watts, an increase of the mean power of 1.0 
watts, and an increase in the input frequency by a factor of two from the 
baseline case. There is heat transfer enhancement for all three heater 
locations. The mean cyclic temperature response is higher than the case just 
discussed, see Figures IA-IC, however, the maximum cyclic temperature 
responses are approximately the same. The amplitude of the temperature 


response is less than for the baseline case. 
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Figures DG-DI present results for mean input of 0.5 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These result from a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.5 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. There is heat transfer enhancement for all three heater locations 
with significant heat transfer enhancement for the top heater location only. 
The amplitude of the temperature response is similar to that of the baseline 
case. 

Figures EG-EI present results for a mean input of 1.0 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These parameters are a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.0 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. There is heat transfer enhancement for all three heater 
locations. The amplitude of the temperature response is similar to that of the 
baseline case. 

Figures GG-GI present results for mean input of 3.0 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These parameters are a 
decrease in the amplitude of 0.5 watts, an increase of the mean power of 1.0 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. There is heat transfer enhancement for all three heater 
locations. The amplitude of the temperature response is similar to that of the 


baseline case. 
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Figures IG-I] presents results for mean input of 3.0 watts, 
amplitude of 2.0 watts and a frequency of 0.025 Hz. These values are an 
increase in the amplitude of 1.0 watts, an increase of the mean power of 1.0 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. There is heat transfer enhancement for all three heater 
locations. The amplitude of the temperature response is much greater than 
that of the baseline case. 

In examining these cases the amplitude to mean input power 
ratio was the dominant factor in producing heat transfer enhancement. 
Increase in the frequency of the input power provided some heat transfer 
enhancement, but more importantly it greatly influenced the amplitude of the 
temperature response. Higher mean power inputs provide some heat transfer 
enhancement. However they were a less dominant factor than the amplitude 


to mean power input ratio. 


C. DISCUSSION OF SQUARE WAVE DATA 

As was the case for the triangular wave data, the base line profile had a 
mean value power input of 2.0 watts, amplitude of 1.0 watts, and a frequency 
of 0.050 Hz, see Figures KD-KF. 

Note the spikes in the power pulsation. These spikes were prevalent 
throughout the square wave data. The exact cause is inconclusive, however, it 


is believed that the spikes are due to the power supply hunting for the proper 
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voltage to match the control signal supplied by the programmer. The control 
signal produced by the programmer is the same as the desired output voltage. 
In examining the raw data, these spikes are caused at isolated points, in some 
cases two points out of a total of nine hundred data points for each run. These 
spikes did not seem to have an impact on the temperature responses. 

The amplitude of the cyclic temperature response for the baseline case is 
14.8% of the mean for the top heater location, 17.3% for the middle heater 
location and 18.4% for the bottom heater location. This is for a power 
amplitude 50% of the mean input. The amplitude of the temperature response 
for the baseline case of the triangular wave power pattern range from 6.7% at 
the top heater location to 16.8 % for the bottom heater location. The frequency 
of the cyclic temperature response matches that of the input power pulsation. 
There is a phase lag in the cyclic temperature response, but not as large as in 
the triangular power pattern. The mean cyclic temperature Pogo coincides 
with the steady temperature response indicating no heat transfer enhancement 


due to power pulsation. 


1. Effects of Amplitude of Power Pulsation 
Figures LD-LF present results for mean input of 2.0 watts, 
amplitude of 2.0 watts and a frequency of 0.050 Hz. This is an increase of 1.0 
watt to the baseline case. There is only a slight improvement in the Heat 
transfer characteristics. The amplitude of the cyclic temperature response is 


increased to 32.6% for the top heater, 42.6% for the middle heater and 55.7% 
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for the bottom heater. The maximum value of the cyclic temperature response 
for the middle heater is actually higher than that of the top heater. 

Figures JD-JF present results for mean input of 2.0 watts, amplitude 
of 0.5 watts and a frequency of 0.050 Hz. This is a decrease in the amplitude 
of 0.5 watts compared to the baseline case. There is no heat transfer 
enhancement for this case. The mean cyclic temperature response coincides 
with the steady temperature response for all three heater locations. The 
amplitude of the cyclic temperature response is decreased significantly to 
8.35% for the top heater location, 10.7% for the middle heater location, and 
10.2% for the bottom heater location. Once again, the maximum cyclic 
temperature response for the middle heater is higher than that of the top 


heater. 


2. Effects of Time Period of Pulsation 

Figures KA-KC present results for mean input of 2.0 watts, 
amplitude of 1.0 watts and a frequency of 0.100 Hz. This is an increase in the 
frequency of pulsation from the baseline case by a factor of two. There is no 
heat transfer enhancement for this case. The amplitude of the temperature 
response for the top heater location drops to 7.8% of the mean, 12.4% for the 
middle heater location, and 10.5% for the bottom heater location. The 
frequency of the cyclic temperature response matches that of the pulsating 
power input. The maximum cyclic temperature response for the middle heater 


location is slightly larger than that of the top heater. 
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Figures KG-KI present results for mean input of 2.0 watts, 
amplitude of 1.0 watts and a frequency of 0.025 Hz. This is a decrease in the 
frequency of pulsation from the baseline case by a factor of two. There is only 
a slight heat transfer enhancement for the top and middle heater locations and 
no enhancement for the bottom heater location. The amplitude of the 
temperature response of the top heater location jumps to 24.5%, the middle 
heater location to 27.5%, and the bottom heater location to 26.3%. The middle 
heater still has a slightly larger maximum cycle temperature than that of the 
top heater. Also noted was a drop in temperature at the end of each cycle. This 


effect is most prominent for the middle heater, Figure KH. 


3. Effects of the Mean Input Power 

Figures QD-QF present results for mean input of 3.0 watts, 
amplitude of 1.0 watts and a frequency of 0.050 Hz. This is an increase of 1.0 
watt to the baseline case. The temperature amplitude for the middle and 
bottom heaters are approximately the same, while the top heater amplitude 
has decreased. The frequency of the cyclic temperature response matches that 
of the power input. The mean cyclic temperature response coincides with the 
steady temperature responses. 

Figures OD-OF present results for mean input of 1.0 watts, 
amplitude of 1.0 watts and a frequency of 0.050 Hz. This is a decrease of 1.0 
watt to the baseline case. By decreasing the mean power input to 1.0 watts, 


the amplitudes of the temperature response for all three heater locations are 
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approximately the same as in the base line case. The frequency response was 
unchanged. The mean of the cyclic temperature response is only slightly lower 
than the temperature response to the steady input for the top and middle 


heaters locations. 
4. Combing of Effects 


a. Varying Amplitude and Frequency 

Figures LA-LC present results for mean input of 2.0 watts, 
amplitude of 2.0 watts and a frequency of 0.100 Hz. These figures are an 
increase in the amplitude of 1.0 watts and an increase of the input frequency 
by a factor of two from the baseline case. There is a slight heat transfer 
enhancement for all three heater locations. The amplitude of the temperature 
response is larger than the baseline case. 

Figures LG-LI present results for mean input of 2.0 watts, 
amplitude of 2.0 watts and a frequency of 0.025 Hz. These figures are an 
increase in the amplitude of 1.0 watts and a decrease of the input frequency 
by a factor of two from the baseline case. For all three heater locations there 
is a slight heat transfer enhancement. Also, there is a large increase in the 
amplitude of the temperature response as compared to the baseline case. 

Figures JA-JC present results for mean input of 2.0 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These parameters are a 


decrease in the amplitude of 0.5 watts and an increase of the input frequency 


4] 


by a factor of two from the baseline case. Only the top heater location 
demonstrates a slight enhancement to the heat transfer. As expected, the 
amplitude of the temperature response is less than the baseline case. 

Figures JG-JI present results for mean input of 2.0 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These values are a 
decrease in the amplitude of 0.5 watts and a decrease of the input frequency 
by a factor of two from the baseline case. No heat transfer enhancement for 
any of the three heater locations is noted. However, there is a dramatic drop 
in temperature at the end of each cycle, as was discussed early. 

No significant enhancement to the heat transfer, as defined 
early, was found for any of the runs performed. The drop in temperature at the 
end of each cycle was apparent for amplitudes of 0.5 watts and an input 
frequency of 0.025 Hz. As for the triangular wave power pattern data, 
frequency of the input power had a influence over the amplitude of the 
temperature response. However, the amplitude to mean power input ratio did 
not bring about the dramatic effects seen with the triangular wave power 


pattern data. 


6b. Varying Amplitude and Mean Input Power 
Figures RD-RF present results for mean input of 3.0 watts, 
amplitude of 2.0 watts and a frequency of 0.050 Hz. These values are an 
increase in the amplitude of 1.0 watts and an increase of the mean power of 


1.0 watts from the baseline case. Only slight heat transfer enhancement can 
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be seen at all three heater locations. The amplitude of the temperature 
response is greater than the baseline case. 

Figures PD-PF present results for mean input of 3.0 watts, 
amplitude of 0.5 watts and a frequency of 0.050 Hz. These values are a 
decrease in the amplitude of 0.5 watts and an increase of the mean power of 
1.0 watts from the baseline case. No heat transfer enhancement was present 
for any of the three heater locations. For the top heater location there is a 
degrade in the heat transfer effects. The amplitude of the temperature 
response is less than the baseline case. 

Figures ND-NF presents results for mean input of 1.0 watts, 
amplitude of 0.5 watts and a frequency of 0.050 Hz. These parameters are a 
decrease in the amplitude of 0.5 watts and a decrease of the mean power of 1.0 
watts from the baseline case. All three heater locations demonstrated a slight 
enhancement to the heat transfer. The amplitude of the temperature response 
is less than the baseline. 

Figures MD-MF present results for mean input of 0.5 watts, 
amplitude of 0.5 watts and a frequency of 0.050 Hz. These parameters are a 
decrease in the amplitude of 1.5 watts and a decrease of the mean power of 0.5 
watts from the baseline case. No heat transfer enhancement for any of the 
three locations is seen. The amplitude of the temperature response is less than 


the baseline. 


43 


The effect of the amplitude of the power input was to increase 
the amplitude of the temperature response. The dramatic effects seen with 


triangular wave data were not seen here. 


c. Varying Frequency and Mean Input Power 

Figures QA-QC present results for mean input of 3.0 watts, 
amplitude of 1.0 watts and a frequency of 0.100 Hz. These values are an 
increase in the mean power input of 1.0 watts and an increase of the input 
frequency by a factor of two from the baseline case. None of the three of the 
heater locations demonstrated any heat transfer enhancement. The amplitude 
of the temperature response is less than the baseline case. 

Figures QG-QI present results for mean input of 3.0 watts, 
amplitude of 1.0 watts and a frequency of 0.025 Hz. These parameters are an 
increase in the mean power input of 1.0 watts and a decrease of the input 
frequency by a factor of two from the baseline case. Only a slight heat transfer 
enhancement is seen for all three heater locations. There was a drop in 
temperature at the end of each cycle. The amplitude of the temperature 
response was greater than the baseline case. 

Figures OA-OC present results for mean input of 1.0 watts, 
amplitude of 1.0 watts and a frequency of 0.100 Hz. These values are a 
decrease in the mean power input of 1.0 watts and an increase of the neat 


frequency by a factor of two from the baseline case. A slight heat transfer 
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enhancement was observed for all three heater locations. The amplitude of the 
temperature response is less than the baseline case. 

Figures OG-OI present results for mean input of 1.0 watts, 
amplitude of 1.0 watts and a frequency of 0.025 Hz. These parameters are a 
decrease in the mean power input of 1.0 watts and a decrease of the input 
frequency by a factor of two from the baseline case. A slight heat transfer 
enhancement is observed at all three heater locations. The amplitude of the 
temperature response is greater than the baseline case. 

The last two cases discussed both had a mean power input of 1.0 
watts and amplitude of 1.0 watts. Both of these cases displayed some heat 
transfer enhancement, implying that the amplitude to mean power input ratio 
of one does have some impact upon the heat transfer characteristics. However, 
this enhancement is small in comparison to the triangular wave data and only 
applies for mean power inputs equal to or greater than 1.0 watts. The large 
temperature drop at the end of the cycle found in Figures OG-OI, aided in the 
heat transfer enhancement. However, the maximum cyclic temperature 


response was not lessened. 


d. Other Variations of Parameters 
Figures MA-MC present results for mean input of 0.5 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These values are a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.5 


watts, and an increase in the input frequency by a factor of two from the 
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baseline case. No heat transfer enhancement is seen for this case. The 
amplitude of the temperature response is less than the baseline case. Also, the 
mean cyclic temperature response is higher than the triangular wave data for 
identical parameters. 

Figures NA-NC present results for mean input of 1.0 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These values are a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.0 
watts, and an increase in the input frequency by a factor of two from the 
baseline case. There was a slight enhancement for all three heater locations. 
The amplitude of the temperature response is less than the baseline case. 

Figures RA-RC present results for mean input of 3.0 watts, 
amplitude of 2.0 watts and a frequency of 0.100 Hz. These values are an 
increase in the amplitude of 1.0 watts, an increase of the mean power of 1.0 
watts, and an increase in the input frequency by a factor of two from the 
baseline case. A slight heat transfer enhancement, approximately the same as 
in the previous case, is seen for all three heater locations. The amplitude of the 
temperature response is less than the baseline case for the top heater location, 
greater for the middle heater location, and approximately the same for the 
bottom heater. Heater configuration C was used for the baseline case while 
heater configuration B is used for this case with different locations of the 


heaters examined. This could explain the difference. 
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Figures PA-PC present results for mean input of 3.0 watts, 
amplitude of 0.5 watts and a frequency of 0.100 Hz. These parameters are a 
decrease in the amplitude of 0.5 watts, an increase of the mean power of 1.0 
watts, and an increase in the input frequency by a factor of two from the 
baseline case. No heat transfer enhancement for any of the three heater 
locations is noted. The amplitude of the temperature response is much less 
than the baseline case. 

Figures MG-MI present results for mean input of 0.5 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These values are a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.5 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. No heat transfer enhancement is found for these conditions. The 
amplitude of the temperature response is approximately the same as the 
baseline case. 

Figures NG-NI present results for a mean input of 1.0 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These values are a 
decrease in the amplitude of 0.5 watts, a decrease of the mean power of 1.0 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. Only slight heat transfer enhancement is seen here. The 
amplitude of the temperature response is approximately the same as in the 
baseline case. The middle heater location is starting to experience a 


temperature drop at the end of the cycle. 
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Figures PG-PI present results for mean input of 3.0 watts, 
amplitude of 0.5 watts and a frequency of 0.025 Hz. These values are a 
decrease in the amplitude of 0.5 watts, an increase of the mean power of 1.0 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. All three heater locations demonstrate a temperature drop at the 
end of each cycle. The middle heater shows the most dramatic drop of the three 
heater locations. There is only a slight heat transfer enhancement seen for 
these three heater locations. The greatest enhancement is seen for the middle 
heater. The amplitude of the temperature response is slightly less than the 
baseline case. 

Figures RG-RI present results for mean input of 3.0 watts, 
amplitude of 2.0 watts and a frequency of 0.025 Hz. These values are an 
increase in the amplitude of 1.0 watts, an increase of the mean power of 1.0 
watts, and a decrease in the input frequency by a factor of two from the 
baseline case. Only a slight heat transfer enhancement are seen for these three 
heater locations. The middle heater shows a small temperature drop at the end 
of the cycle. The amplitude of the temperature response is much greater than 
in the baseline case. 

These experiments reveal a much less pronounced trend for heat 
transfer enhancement then observed for the triangular wave. Large ratios of 
amplitude to mean power level have a minor but positive impact on the heat 


transfer characteristics. Frequency of the power pulsation has an effect on the 
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amplitude of the temperature response, but once again much less than for the 
triangular wave data. Higher mean power inputs influence the large 
temperature drops at the end of each cycle. However, this enhancement does 


not lead to significant heat transfer enhancement over the entire cycle. 


D. NON DIMENSIONAL HEAT TRANSFER CHARACTERISTICS 
Figure 8 presents a comparison of the steady data for the bottom heater 
location with previous investigations. The non dimensional data are calculated 


by defining the Nusselt number, Nu, and Grashof number, Gr as: 


“us Qk 
PAG =a») 


where Q..... was determined by the following equation 


conv 


Q conv = Q 7 Q sond 


Q.ong Was estimated from [Ref. 16] 


Cae gBQL* 
kAv’ 


Along with the present measurements also contained within Figure 8 is the 
applicable steady data from Haukenes [Ref. 5] and Akdeniz [Ref. 17], present 
measurements were obtained with the same test surface as used by there 


investigations. Also seen in Figure 8 is an empirical correlation from [Ref. 15]: 


Nu VAI Gr) 
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Figure 8. Comparison of Steady Data 


The present data are in good agreement with those of Haukenes and Akdeniz 
for the same heater spacing. They are about 20% below the empirical equation. 
The lower Nu’ values are expected due to the reduced heater spacing in the 
present measurements compared to [Ref. 15] which results in larger interaction 
between heaters by substrate conduction. 

Figures 9 and 10 are variations of the average Nusselt number with 
average Grashof number with amplitude of pulsation as a parameter. From the 


cyclic responses of heater temperature and power inputs, mean values were 


50 


Nu 





Figure 9. Average Nusselt Number Versus Average Grashof 
Number for the Top Heater Location, Frequency of 
0.050 Hz for Triangular Wave Power Pattern. 


determined. The average Nusselt number and average Grashof number were 
then determined as previously discussed based on these cycle averaged values. 
Figure 9 is for triangular wave power pattern and Figure 10 is for square wave 
power pattern. All points are for the top heater location and frequency of 
pulsation of 0.050 Hz, and varying mean input power levels. The upper three 
points for Figure 9 are for an amplitude to mean power input ratio of one. This 


dramatically illustrates the effect to the heat transfer characteristics as the 


oO] 


Nu 





100 1000 =a 10000 100000 


Eq Amplitude of 0.5 W >< Amplitudeof10W 4 Amplitude of 2.0W 





Figure 10. Average Nusselt Number Versus Average Grashof 
Number for the Top Heater Location, Frequency of 
0.050 Hz for Square Wave Power Pattern. 


amplitude to mean power input approaches one. There is the expected trend 
of increasing Nusselt number for higher Grashof number which indicates that 
more vigorous flows lead to higher heat transfer rates. 
Unlike the triangular wave power pattern, Figure 9, there is no dramatic 
heat transfer enhancement seen for the square wave power pattern. 
Uncertainty analysis is contained within Appendix E for the above 


calculations. No attempt was made to include the response time of the 


o2 


thermocouple. Also, the uncertainty of the measurements of the physical 
dimensions of the heater, surface area and perimeter, were considered 
negligible compared to the uncertainty of the temperature and power 


measurements. 
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IV. CONCLUSIONS 


In comparing the base line profiles of the two data sets, the square wave 
power pulsation produces less thermal damping effects than the triangular 
wave power pulsation. In addition, the degree of thermal damping effects are 
approximately the same for all the three heaters for the square wave power 
pulsation. Also the phase lag in the cyclic temperature response is more 
pronounced for the triangular wave power pattern. 

With the triangular wave power pattern, significant enhancement to the 
heat transfer was achieved. This was most pronounced when the ratio of 
amplitude to mean power input was the greatest, for the upper heater, for 
increased frequency and for higher mean power inputs. 

The square wave power input was only marginally effective in lowering 
the mean cyclic temperature response below the temperature response for the 
steady mean power level. However, dramatic temperature drops over part of 
the temperature response cycle were found for a frequency of 0.025 Hz, smaller 
amplitudes and higher mean power input values. This drop occurred at the end 
of each cycle. The maximum cyclic temperature response was never below the 


steady temperature response as in the triangular wave pattern. 


o4 


With the ultimate objective of improving heat transfer characteristics for 
electronic components, triangular wave power pattern holds much promise. 
Square wave power patterns only produced marginal results. Specifically, the 
triangular wave pattern power pulsation provided improvements to the heat 


transfer characteristics; 


¢ when the amplitude to mean power input approached one 
¢ at higher mean power inputs (3.0 watts) 


e at higher frequencies (0.100 Hz) 


D0 


V. RECOMMENDATIONS 
The following recommendations are made for further study: 
Vary the pattern of the power input, starting with a sine wave 


pattern. 


Investigate larger frequencies and higher mean power input levels for 
the triangular wave power pattern. 


For square wave power pattern, further investigation is needed to 
determine the reason behind the large drop in temperature at the end 
of the cycle for higher mean input power levels, lower frequencies and 
lower amplitudes. 
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APPENDIX A: LIST OF RUNS 


Number Value Pattern 
aa [2 | 05 | 0100 | Top | Triangular 
[ap [ 2 | 05 | 0.100 | Middle | ‘Triangular 
[ap [2 | 05 | 0.050 | Top | Triangular 
calf 2 | 08 | 0025 | Bottom | ‘Triangular _ 
[pa] 21 ~~ 0.100 ‘Top | ‘Triangular 
pp | ooo ate ingle 
ac | 2 | 1 | 0.100 | Bottom | Triangular _ 
[sp] 2 | 1 | 0050 | Top | Triangular 
se | 2 | 1 | 0050 | Miaale_[ Triangular 
er [2 | 1 | 0.050 [Bottom | Triangular 
[pe] 2 | 1 | 0025 | ‘Top | Triangular 
pa 2 | 1 | 0025 | “Miade | Triangular 
—s_|_2 | 1 | 0025 | Bottom | Triangular 
[ca [2 | 2 | 0100 | Top | ‘Triangular 
[ca | 2 | 2 | 0100 | Middle | Triangular _ 
ec [2 | 2 | 0100 | Bottom | Triangular _ 
[ep] 2 | 2 | 0050 | Top | ‘Triangular 
ce [2 | 2 | 0050 | Middle _| Triangular 
—cr_ [| 2 | 2 | 0.050 | Bottom | Triangular 
[cc | 2 | 2 | 0025 | Top | ‘Triangular 
[cH [ 2 | 2 | 0.035 | Miaale_[ “Triangular 
af 2 | 2 | 0.025 | Bottom | ‘Triangular 
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aoe ‘| Amplitude | Frequency Heater Power 
Number Value Pattern 
a Pen 
a a 
[pp [05 | 05 | 0050 | Top | ‘Triangular 
[pF [_05 | 05 | 0.050 | Botwom | Triangular 
A 
or [08 
Cet f 1 | 05 | 0025 | Bottom | Triangular 
[Fa [aad n00 |Top | Triana 
[ref] 1 | 0.100 | Midale [Triangular 
[re [0.100 | Bottom | Triangular 
[wp] 1) 0050 |Top | Triangular 
A 
Cre [ 1 | _0.050_ | Roteom | Triangular 
[re [YY 0025 |Top | Triangular 
[ref 0025 [Maat [Triangular 
Tt ft 0.025 Bottom | Triangular 






















































og 


Number Value Pattern 
ca os 
ac 

ap [3 [05 | 0.050 | Top | Tiangular 

[ce [3 | 05 | 0.050 | Middle [Triangular 

[_er_| 3 | 05 | 0.050 | Bottom | Triangular 

[ee [3] 05 | 0025 | Top | Triangular 

[en [3 [05 | 0025 | Miaate | Triangular 

[_ar_| 3 [05 | 0.025 | Bottom | Triangular _ 

[Ha [31 | 0100 | Top | Triangular 

[ap [3 [1 | 0100 | Middle [Triangular _ 

_ac__[ 3 | 1 | 0.100 | Bottom | Triangular 

[ap [a 0050 |Top | Triangolar 

HE [10.050 | idate [Triangular 

[0.050 | Bottom | Triangular 

[ug [3 | 1 | 0025 | Top | Triangular _ 

[auf 3a 0025 | Migate [Triangular | 

Car | 3 | 1 | __0.025 | Bottom | Triangular 

[ta [0100 |Top |) Triangular | 

[ps ot00 | Migate | Triangular | 

1D 


G2 
ee 
ee ees eee 


CO} Co 
bo | bo 


— 
eo) 


= 


0.025 | Mid 
0.028 


= 
e9) 


—_— 


Number Value Pattern 
a a 
at] 2 | 05 | 0.025 [Bottom | Square _ 
A 
[Kp [| 0300 | Middle Square 
—Kc__| 2 | 1 | 0100] Bottom | Square 
| kD | 2 | 2 | 0050 | Top | Square 
| KE {| 2 [{ 1 {| 0.050 | Middle {| Square 
KF {| 2 {1 t_ 0.050 | _Bottom | Square 
CT iat 
[Ke] 1 | 0.025 | Maat | Sauare 
«| 2 | 1 | 0025 | “Bottom | Square _ 
a 
[ip [2 | 000 | Middle | Square 
ic [2 | 2 | 01100 | Bottom | Square _ 
up| oso [Top | Square 
[te [2 [2 | 0.050 | Middle | Square 
aur [ 2 | 2 | 0.050 | “Bowtom [Square _ 
7 
[tH [2 [3 [0.025 | “Miaate [Square 
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Mean Frequency Power 
Number Value Pattern 
0.5 05 | 0300 | ee 
| MB | 05 =| 05 | 0.100 | Middle | Square 
| ME | 95 | 05 |__ 0-100 | Bottom | Square 
[wo [05 | 0s [0050 [Top 
[np [1 * os | 0.100 | Middle 
a ee 
a 
oc fa] | 0.200 | Bottom 
KS 
[oz | tf 0.050 [Maat 
orf 1] 1 __| 0.050 [Bottom 
A 
[on [a 0.085 ate 
ort 0.025 [Bottom 



































Square 







Square 
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Figure | Amplitude Frequency Heater Power 
Number Pattern 


[pp [3 [05 | 0.100 | ia 
0.050 | Middle 
0.025 
0.025 | Middle 
Middl 
a 
0.050 | Middle 

ac 0.028 

a Middle 
[eps | 0100 | miaate 
0.050 | Middle 




















































0.035 | Midd 
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APPENDIX B: REDUCE DATA 
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APPENDIX C: 
TEMPERATURE AND POWER ACQUISITION PROGRAM 
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'!! TEMPERATURE AND POWER TRANSIENT RESPOSE PROGRAM ae 


TUOPPPOPTPPP PTET TT OPP ETOP TT ee PT eee eee eee terre rer eeroee ere eeeey 


eeeeeeeoee Cee eee esteoeeeeee eee esteeeesee ese seston eeerteseeeeeeeteeenerece 


OPTION BASE 1 

REAL T(2000),Volts(2000),Pow(2000),V(2000),Tyme(2000),Tsec(2000) 
REAL Vminus(2000), Vplus(2000),Tminus(2000),Tplus(2000) 
REAL V46(2000),V47(2000),V48(2000),V49(2000),V50(2000),V51(2000) 
REAL Temp(3), Vamb(3) 

FOR K=1 TO 2000 

HGS=0. 

Volts(K)=0. 

Pow(K)=0. 

V(K)=0. 

Tyme(K)=0. 

Tsec(K)=0. 

Vminus(K)=0. 

Vplus(K)=0. 

Tminus(K)=0. 

Tplus(K)=0. 

V46(K)=0. 

V47(K)=0. 

V48(K)=0. 

V49(K)=0. 

V50(K)=0. 

VolhKj=0. 

NEXT K 

' 


PRINTER IS CRT 

! 

INPUT "PLEASE INPUT THE NUMBER OF THE HEATER OF INTEREST".N 
t 


IF N=1 THEN 


PRINT "THIS THERMOCOUPLE IS BAD, HOWEVER I CAN GIVE YOU THE" 
PRINT "TEMPERATURE FOR HEATER NUMBER 31 WHICH 


GEOMETRICALLY" 


350 
360 
370 
380 
390 
400 
410 
420 


PRINT "SYMETERICAL TO HEATER 1° 
PRINT * " 
PRINT "PRESS Continue TO RESUME" 
PAUSE 

END IF 


! 


le N=3 THEN 


PRINT "THIS THERMOCOUPLE IS BAD HOWEVER I CAN GIVE YOU THE" 


2a 


IS 


430 
440 
450 
460 
470 
480 
490 
300 
210 
020 
030 
540 
5900 
5960 
O70 
5980 
090 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
120 
730 
740 
750 
760 
770 
780 
720 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 


PRINT "AVERAGE TEMPERATURE OF THE THERMOCOUPLES ABOVE AND" 
PRINT "BELOW" 
PRINT " " 
PRINT "PRESS Continue TO RESUME" 
PAUSE 
END IF 
IF N=6 THEN 
PRINT "THIS THERMOCOUPLE IS BAD HOWEVER I CAN GIVE YOU THE" 
PRINT “AVERAGE TEMPERATURE OF THE THERMOCOUPLES ABOVE AND" 
PRINT "BELOW" 
PRINT °° 
PRINT "PRESS Continue TO RESUME" 
PAUSE 
END IF 


IF N=22 THEN 
PRINT "THIS THERMOCOUPLE IS BAD HOWEVER I CAN GIVE YOU THE" 
PRINT "AVERAGE TEMPERATURE OF THE THERMOCOUPLES ABOVE AND" 
PRINT "BELOW" 
PRENG 
PRINT "PRESS Continue TO RESUME" 
PAUSE 

END IF 


IF N=29 THEN 
PRINT "THIS THERMOCOUPLE IS BAD HOWEVER I CAN GIVE YOU THE" 
PRINT “AVERAGE TEMPERATURE OF THE THERMOCOUPLES ABOVE AND" 
PRINT "BELOW" 
PRINT " " 
PRINT "PRESS Continue TO RESUME" 
PAUSE 
END IF 
! 
IF N=37 THEN 
PRINT "THIS THERMOCOUPLE IS BAD HOWEVER I CAN GIVE YOU THE" 
PRINT "AVERAGE TEMPERATURE OF THE THERMOCOUPLES ABOVE AND" 
PRINT "BELOW" 
PRINT " " 
PRINT "PRESS Continue TO RESUME" 
PAUSE 
END IF 
INPUT "PLEASE INUPT THE CYCLE TIME IN SECONDS',Np 
INPUT "PLEASE INUPT THE MAX AND MIN VOLTAGE",Maxv,Minv 
INPUT "PLEASE INUPT THE OUTPUT OF THE SIDE COLUMNS IN WATTS',Scp 
IF Np<10 THEN GOTO 860 
Wx=(Np/50)-1.0 
IF Wx<=0 THEN Wx=0. 
M=INT((60/Np)*250)+500 
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930 

940 

950 

960 

970 

980 

990 

1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
T1Z0 
1130 
1140 
1150 
1160 
P70 
1180 
1190 
1200 
1210 
1220) 
1230 
1240 
1250 
1260 
1278 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 


IF M>900 THEN M=900 


! 


Ti=TIMEDATE 


FOR I=1 TOM 


Mu tPUT 709: RST 


IF N=) THEN 


OUTPUT 709:"CONFMEAS DCV,210,USE 0" 
ENTER 709;V(1) 

OUTPUT 709;"CONFMEAS DCV,315,USE 0" 
ENTER 709;Volts(I) 

OUTPUT 709;'CONFMEAS DCV,316,USE 0" 
ENTER 709;V46(1) 

GOTO 5590 


END IF 


te N=2 THEN 


OUTPUT 709;"CONFMEAS DCV,101,USE 0" 
ENTER 709;V(1D) 

OUTPUT 709;"CONFMEAS DCV,314,USE 0" 
ENTER 709;Volts(1) 

OUTPUT 709;"CONFMEAS DCV,316,USE 0° 
ENTER 709;V46(1) 

GOTO 5590 


END IF 


IF N=3 THEN 


OUTPUT 709;"CONFMEAS DCV,101,USE 0° 
ENTER 709;Vminus(I) 

OUTPUT 709;'CONFMEAS DCV,103,USE 0° 
ENTER 709;Vplus(1) 

OUTPUT 709;"CONFMEAS DCV,313,USE 0" 
ENTER 709;Volts(D 

OUTPUT 709;"CONFMEAS DCV,316,USE 0° 
ENTER 709;V46(1) 

GOTO 5590 


END IF 


IF N=4 THEN 


OUTPUT 709;"CONFMEAS DCV,103,USE 0° 
ENTER 709;V(I) 

OUTPUT 709;'CONFMEAS DCV,312,U0SE 0° 
ENTER 709;Volts(D) 

OUTPUT 709;"CONFMEAS DCV,316,USE 0° 
ENTER 709;V46(1) 

GOTO 5590 


END IF 


IF N=5 THEN 
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1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 


OUTPUT 709;"CONFMEAS DCV,104,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,311,USE 0" 
ENTER 709;Volts(D 
OUTPUT 709;"CONFMEAS DCV,316,USE 0" 
ENTER 709;V46(1) 
GOTO 5590 

END IF 

IF N=6 THEN 
OUTPUT 709;"CONFMEAS DCV,104,USE 0" 
ENTER 709;Vminus(I) 
OUTPUT 709;"CONFMEAS DCV,106,USE 0" 
ENTER 709;Vplus(I) 
OUTPUT 709;"CONFMEAS DCV,310,USE 0" 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,316,USE 0" 
ENTER 709;V46(1) 
GOTO 5590 

END IF 

' 

IE Ne? THEN 
OUTPUT 709;"CONFMEAS DCV,106,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,309,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,316,USE 0" 
ENTER 709;V46(1) 
GOTO 5590 

END IF 

IF N=8 THEN 
OUTPUT 709;"CONFMEAS DCV,100,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,308,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;'CONFMEAS DCV,316,USE 0" 
ENTER 709;V46(1) 
GOTO 5590 

END IF 

' 

IF N=9 THEN 
OUTPUT 709;"CONFMEAS DCV,100,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,403,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,400,USE 0" 
ENTER 709;V47(1) 
GOTO 5590 

END IF 
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1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
ZO 
2130 
2140 
ZO 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 


IF N=10 THEN 
OUTPUT 709;"CONFMEAS DCV,109,USE 0° 
ENTER 709;V(1) 
OUTPUT 709;'CONFMEAS DCV,404,USE 0" 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,400,USE 0" 
ENTER 709;V47(1) 
GOTO 5590 

END IF 


IF N=11 THEN 
OUTPUT 709:'CONFMEAS DCV,110,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;'CONFMEAS DCV,405,USE 0" 
ENTER 709;Volts(I) 
OUTPUT 709;'CONFMEAS DCV,400,USE 0° 
ENTER 709;V47(1) 
GOTO 5590 

END IF 

IF N=12 THEN 
OUTPUT 709;'CONFMEAS DCV,111,USE 0" 
ENTER 709;V(1) 
OUTPUT 709:"CONFMEAS DCV,406,USE 0" 
ENTER 709:Volts(1) 
OUTPUT 709;"CONFMEAS DCV,400,USE 0° 
ENTER 709;V47(1) 
GOTO 5590 

END IF 

] 

IF N=13 THEN 
OUTPUT 709;"CONFMEAS DCV,112,USE 0° 
ENTER 709:V(1) 
OUTPUT 709;"CONFMEAS DCV,407,USE 0" 
ENTER 709-;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,400,USE 0° 
ENTER 709:V47(1) 
GOTO 5590 

END IF 

IF N=14 THEN 
OUTPUT 709:"CONFMEAS DCV,113,USE 0° 
ENTER 709:V(1) 
OUTPUT 709;'CONFMEAS DCV,408,USE 0° 
ENTER 709-Volts(I) 
OUTPUT 709;"CONFMEAS DCV,400,USE 0° 
ENTER 709;V47(1) 
GOTO 5590 

END IF 
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2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2930 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2916 
2920 


IF N=15 THEN 
OUTPUT 709;"CONFMEAS DCV,114,USE 0" 
ENTER 709;V(I) 
OUTPUT 709;"CONFMEAS DCV,409,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,400,USE 0” 
ENTER 709;V47(1) 
GOTO 5590 

END IF 

IF N=16 THEN 
OUTPUT 709;"CONFMEAS DCV,115,USE 0” 
ENTER 709; V(I) 
OUTPUT 709;"CONFMEAS DCV,410,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,401,USE 0" 
ENTER 709;V48(I) 
GOTO 5590 

END IF 

IF N=17 THEN 
OUTPUT 709:"CONFMEAS DCV,116,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,411,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,401,USE 0" 
ENTER 709;V48(1) 
GOTO 5590 

END IF 

' 


IF N=18 THEN 
OUTPUT 709;"CONFMEAS DCV,117,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,412,USE 0" 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,401,USE 0" 
ENTER 709;V48(1) 
GOTO 5590 

END IF 

' 

IF N=19 THEN 
OUTPUT 709;"CONFMEAS DCV,118,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,413,USE 0" 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,401,USE 0" 
ENTER 709;V48(I) 
GOTO 5590 

END IF 
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2930 
2940 
2950 
2960 
2070 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 


! 


VEIN =20 THEN 
OUTPUT 709;"CONFMEAS DCV,119,USE 0" 
ENTER 709;V(1) 
OUTPUT 709:"CONFMEAS DCV,414,USE 0" 
ENTER 709:Volts(I) 
OUTPUT 709;"'CONFMEAS DCV,401,USE 0" 
ENTER 709;V48(1) 
GOTO 5590 

END IF 

IF N=21 THEN 
OUTPUT 709;"CONFMEAS DCV,200,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,415,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709:"CONFMEAS DCV,401,USE 0" 
ENTER 709:V48(1) 
GOTO 5590 

END IF 

IF N=22 THEN 
OUTPUT 709:"CONFMEAS DCV,200,USE 0" 
ENTER 709:Vminus(1) 
OUTPUT 709;"CONFMEAS DCV,202,USE 0" 
ENTER 709:;VplusiI) 
OUTPUT 709:"CONFMEAS DCV,416,USE 0° 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,401,USE 0° 
ENTER 709;V48(1) 
GOTO 5590 

END IF 

IF N=23 THEN 
OUTPUT 709;"CONFMEAS DCV,202,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,417,USE 0° 
ENTER 709:Volts(1) 
OUTPUT 709."CONFMEAS DCV,401,USE 0° 
ENTER 709-;V48(]) 
GOTO 5590 

END IF 


IF N=24 THEN 
OUTPUT 709;"CONFMEAS DCV,203,USE 0° 
ENTER 709;V(1) 
OUTPUT 709;'CONFMEAS DCV,503,USE 0° 
ENTER 709;Volts(1) 
OUTPUT 709;"'CONFMEAS DCV,500,USE 0° 
ENTER 709;V49(1) 
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3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3010 
3520 
3030 
3540 
35950 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 


GOTO 5590 

END IF 

! 

IF N=25 THEN 
OUTPUT 709;"CONFMEAS DCV,204,USE 0" 
ENTER 709;V‘1) 
OUTPUT 709;"CONFMEAS DCV,504,USE 0" 
ENTER 709:Volts(I) 
OUTPUT 709:"CONFMEAS DCV,500,USE 0" 
ENTER 709:V49(1) 
GOTO 5590 

END IF 

t 


IF N=26 THEN 
OUTPUT 709;"CONFMEAS DCV,205,USE 0° 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,505,USE 0" 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,500,USE 0" 
ENTER 709;V49(1) 
GOTO 5590 

END IF 

t 

IF N=27 THEN 
OUTPUT 709;"CONFMEAS DCV,206,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,506,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,500,USE 0" 
ENTER 709;V49(1) 
GOTO 5590 

END IF 


IF N=28 THEN 
OUTPUT 709;"CONFMEAS DCV,207,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,507,USE 0” 
ENTER 709;Volts(J) 
OUTPUT 709;"CONFMEAS DCV,500,USE 0° 
ENTER 709;V49(T) 
GOTO 5590 

END IF 

IF N=29 THEN 
OUTPUT 709;"CONFMEAS DCV,207,USE 0" 
ENTER 709;Vminus(I) 
OUTPUT 709;"CONFMEAS DCV,209,USE 0" 
ENTER 709;Vplus(1) 
OUTPUT 709;"CONFMEAS DCV,508,USE 0" 
ENTER 709;Volts(1) 
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3930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 


OUTPUT 709;"CONFMEAS DCV,500,USE 0" 
ENTER 709:V49(I) 
GOTO 5590 

END IF 

i 

IF N=30 THEN 
OUTPUT 709;"CONFMEAS DCV,209,USE 0" 
ENTER 709;V(I) 
OUTPUT 709;"CONFMEAS DCV,509,USE 0" 
ENTER 709:Volts(1) 
OUTPUT 709;"CONFMEAS DCV,500,USE 0" 
ENTER 709;V49(1) 
GOTO 5590 

END IF 


im N=3i THEN 
OUTPUT 709;"CONFMEAS DCV,210,USE 0” 
ENTER 709:V(1) 
OUTPUT 709;"CONFMEAS DCV,510,USE 0" 
ENTER 709:Volts(D 
OUTPUT 709;"CONFMEAS DCV,501,USE 0" 
ENTER 709;V50(1) 
GOTO 5590 

END IF 

PpoN=o2 THEN 
OUTPUT 709;"CONFMEAS DCV,211,USE 0" 
ENTER 709:V(1) 
OUTPUT 709:"CONFMEAS DCV,511,USE 0" 
ENTER 709:Volts(1) 
OUTPUT 709;"CONFMEAS DCV,501,USE 0" 
ENTER 709;V50(1) 
GOTO 5590 

END IF 

He N=33 THEN 
OUTPUT 709;"CONFMEAS DCV,212,USE 0° 
ENTER 709:ViJ) 
OUTPUT 709;"CONFMEAS DCV,512,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,501,USE 0" 
ENTER 709;V50(1) 
GOTO 5590 

END IF 

IF N=34 THEN 
OUTPUT 709: CONEMBEAS DOV .213,UsE 0. 
ENTER 709;V(1) 
OUTPUT 709:"CONFMEAS DCV,513,USE 0° 
ENTER 709:Volts(1) 
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4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 


OUTPUT 709;"CONFMEAS DCV,501,USE 0° 
ENTER 709;V50(1) 
GOTO 5590 

END IF 

IF N=35 THEN 
OUTPUT 709;"CONFMEAS DCV,214,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;'CONFMEAS DCV,514,USE 0° 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,501,USE 0° 
ENTER 709;V50(1) 
GOTO 5590 

END IF 

IF N=36 THEN 
OUTPUT 709;"CONFMEAS DCV,215,USE 0° 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,515,USE 0° 
ENTER 709;Volts(1) 
OUTPUT 709;°CONFMEAS DCV,501,USE 0" 
ENTER 709;V50(1) 
GOTO 5590 

END IF 


IF N=37 THEN 
OUTPUT 709;"CONFMEAS DCV,215,USE 0" 
ENTER 709;Vminus(J) 
OUTPUT 709;"'CONFMEAS DCV,217,USE 0° 
ENTER 709;Vplus(J) 
OUTPUT 709;"CONFMEAS DCV,516,USE 0° 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,501,USE 0° 
ENTER 709;V50(1) 
GOTO 5590 

END IF 


IF N=38 THEN 
OUTPUT 709;"CONFMEAS DCV,217,USE 0° 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,517,USE 0° 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,501,USE 0° 
ENTER 709;V50(1) 
GOTO 5590 

END IF 

IF N=39 THEN 
OUTPUT 709;"CONFMEAS DCV,218,USE 0° 
ENTER 709;V(1) 
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4930 
4940 
4950 
4960 
4970 
4980 
4990 
2000 
9010 
9020 
5030 
5040 
9050 
5060 
5070 
5080 
9090 
3100 
9110 
9120 
59130 
5140 
0150 
9160 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 
5260 
5270 
5280 
3290 
5300 
50310 
9320 
5330 
5340 
5350 
5360 
0370 
5380 
5390 
5400 
5410 
5420 


OUTPUT 709;"CONFMEAS DCV,602,USE 0" 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,600,USE 0" 
RENTER 709-Vold 
GOTO 5590 

END IF 


IF N=40 THEN 
OUTPUT 709;"CONFMEAS DCV,219,USE 0" 
ENTER 709:V(1) 
OUTPUT 709;"CONFMEAS DCV,603,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,600,USE 0" 
ENTER 709;V51(1) 
GOTO 5590 

END IF 


IF N=41 THEN 
OUTPUT 709;"CONFMEAS DCV,300,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,604,USE 0” 
ENTER 709;Volts(I) 
OUTPUT 709;"CONFMEAS DCV,600,USE 0” 
ENTER 709;V51(I) 
GOTO 5590 

END IF 

] 


IF N=42 THEN 
OUTPUT 709:"CONFMEAS DCV,301,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,605,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,600,USE 0" 
ENTER 709;V51(1) 
GOTO 5590 

END IF 


IF N=43 THEN 
OUTPUT 709;"CONFMEAS DCV,302,USE 0" 
ENTER 709;V(1) 
OUTPUT 709;"CONFMEAS DCV,606,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,600,USE 0" 
ENTER 709;V51(1) 
GOTO 5590 

END IF 

IF N=44 THEN 
OUTPUT 709:"CONFMEAS DCV,303,USE 0" 
ENTER 709:V(1) 
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5430 
5440 
5450 
5460 
5470 
5480 
5490 
9000 
9010 
9020 
9530 
5040 
5550 
0560 
0070 


5580 ! 


5590 
5600 
5610 
5620 
5630 
5640 
5650 
5660 
9670 
9680 
59690 
5700 
5710 
0720 
0730 
5740 
5750 
5760 
5770 
5780 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
0860 
5870 
5880 
0890 
5900 
59910 
0920 


OUTPUT 709;'CONFMEAS DCV,607,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"CONFMEAS DCV,600,USE 0" 
ENTER 709;V51(1) 
GOTO 5590 

END IF 

IF N=45 THEN 
OUTPUT 709;"CONFMEAS DCV,304,USE 0" 
ENTER 709;V(I) 
OUTPUT 709;"CONFMEAS DCV,608,USE 0" 
ENTER 709;Volts(1) 
OUTPUT 709;"'CONFMEAS DCV,600,USE 0" 
ENTER 709;V51(1) 

END IF 


WAIT Wx 


NEXT I 


! 


Tf=TIMEDATE 

! 

CLEAR SCREEN 
' 


BEEP 81.38,.5 

BEEP 2115.88,.5 
BEEP 3499.34,.5 
BEEP 5208.32,.5 
BEEP 3499.34,.5 
BEEP 2115.88,.5 
BEEP 81.38,.5 


OUTPUT 709;"RST" 

A=.0006797 

B=25825. 1328 

C=607789.2467 

D=21952034.3364 

E=8370810996. 1874 

OUTPUT 709;"CONFMEAS DCV,317-319,USE 0" 
FOR I=1 TO 3 

ENTER 709;Vamb(]) 
Temp(I)=A+(B*Vamb(I))-(C*(Vamb(I)*2))-(D*(Vamb(1)%3))+(E*(Vamb(1)%4)) 
NEXT I 

' 


Tarnbe Genin! 1)+Temp(2)+Temp(3))/3 
] 


OUTPUT 709:"RST" 
FOR I=1 TO M 
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5930 ! 

5940 IF N=3 THEN GOTO 6040 

5950 IF N=6 THEN GOTO 6040 

5960 IF N=22 THEN GOTO 6040 

5970 IF N=29 THEN GOTO 6040 

5980 IF N=37 THEN GOTO 6040 

5990 ! 

6000 T(I=A+(B*V(I))-(C*(V(1)%2))-(D*(V(1)%3))+CE*(V(1)%4)) 

6010 ! 

6020 GOTO 6100 

6030 ! 

6040 Tminus(I)=A+(B* Vminus(I))-(C*(Vminus(I)*2))-(D*(Vminus(1)*3))+(E*(Vminus(1)%4)) 
6050 ! 

6060 Tplus(I)=A+(B*Vplus(1))-(C*(Vplus(I)*2))-(D*(Vplus(1)%3))+(E*(Vplus(1)%4)) 
6070 T(1=(Tminus(1)+Tplus(I))/2 

6080 ! 

6090 ! 

6100 Resist=2.0 

6110 R2=Resist*Resist 

6120 IF N=2 THEN Pow(D=(V46(I)-Volts(1))*(Volts(I)*(Resist+.4)-V46(1)*.4/R2 
6130 IF N=5 THEN Pow(])=(V46(1)-Volts(I))*(Volts(I)*(Resist+.2)-V46(1)*.2)/R2 
6140 IF N=9 THEN Pow(]=(V47(1)-Volts(1))*(Volts(I)*(Resist+.3)-V47(1)*.3)/R2 
6150 IF N=10 THEN Pow(D=(V47(1)-Volts(1))*(Volts(D)*(Resist+.45)-V47(1)*.45)/R2 
6160 IF N=18 THEN Pow()=(V48(1)-Volts(1))*(Volts(1)*(Resist+.55)-V48(1)*.55)/R2 
6170 IF N=23 THEN Pow(1)=(V48(I)-Volts(1))*(Volts(I)*(Resist+.55)-V48(1)*.55)/R2 
6180 IF N=24 THEN Pow(]=(V49(1)-Volts(1))*(Volts(1)*(Resist+. 10)-V49(1)*.10)/R2 
6190 IF N=25 THEN Pow(1D=(V49(1)-Volts(1))*(Volts(1)*(Resist+.30)-V49(1)*.30)/R2 
6200 IF N=32 THEN Pow(I)=(V50(1)-Volts(1))*(Volts(1)*(Resist+.60)-V50(1)*.60)/R2 
6210 IF N=34 THEN Pow(D=(V50(1D-Volts(1))*(Volts(I)*(Resist+.30)- V50(1)*.30)/R2 
6220 IF N=35 THEN Pow(D=(V50(I)-Volts(1))*(Volts(I)*(Resist+.60)-V50(1)*.60)/R2 
6230 IF N=36 THEN Pow(])=(V50(1)-Volts(1))*(Volts()*(Resist+.60)-V50(1)*.60)/R2 
6240 IF N=37 THEN Pow(I)=(V50(1)-Volts(1))*(Volts()*(Resist+.60)-V50(1)*.60)/R2 
6250 IF N=39 THEN Pow(I=(V51(1)-Volts(1I))*(Volts(I)*(Resist+.30)-V51(1)*.30)/R2 
6260 IF N=40 THEN Pow(I)=(V51(1)-Volts(1))*(Volts(1)*(Resist+.30)-V51(1)*.30)/R2 
6270 IF N=42 THEN Pow(I=(V51(1)-Volts(1))*(Volts(1)*(Resist+.55)-V51(1)*.55)/R2 
6280 IF N=43 THEN Pow(D=(V51(1)-Volts(1))*(Volts(I)*(Resist+.45)-V51(1)*.45)/R2 
6290 IF N=44 THEN Pow(I)=(V51(1)-Volts(I))*(Volts(I)*(Resist+.55)-V51(1)*.55)/R2 
6300 IF N=25 THEN Pow(I)=(V49(1)-Volts(1))*(Volts(1)*(Resist+.3)-V49(1)*.3)/R2 
6310 IF N=29 THEN Pow(I=(V49(1)-Volts(1))*(V olts(1)*(Resist+.3)-V49(1)*.3)/R2 
6320 IF N=30 THEN Pow(ID=(V49(1)-Volts(1))*(Volts(I)*(Resist+.3)-V49(1)*.3)/R2 
6330 ! 

6340 Resist=1.972 

6350 R2=Resist* Resist 

6360 ! 

6370 IF N=15 THEN Pow(])=(V47(1)-Volts(1))*(Volts()*(Resist+.2)-V47(I)*.2)/R2 
6380 ! 

6390 Resist=1.9184 

6400 R2=Resist* Resist 

6410 ! 

6420 IF N=14 THEN Pow(D=(V47(I)-Volts(1))*(Volts(1)*(Resist+. 1)-V47(I)*. 1)/R2 
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6430 ! 

6440 Resist=1.886 

6450 R2=Resist*Resist 

6460 ! 

6470 IF N=4 THEN Pow(I)=(V46(1)-Volts(I))*(Volts(1)*(Resist+.1)-V46(1)*.1)/R2 
6480 IF N=6 THEN Pow(I)=(V46(I)-Volts(1))*(Volts(1)*(Resist+.1)-V46(I)*.1)/R2 
6490 ! 

6500 Resist=1.862 

6510 R2=Resist*Resist 

6520 ! 

6530 IF N=27 THEN Pow(1)=(V49(I)-Volts(1))*(Volts(1)*(Resist+.5)-V49(1I)*.5)/R2 
6540 IF N=28 THEN Pow(])=(V49(I)-Volts(1))*(Volts(1)*(Resist+0.)-V49(1)*0.)/R2 
6550 IF N=31 THEN Pow(1)=(V50(1)-Volts(1))*(Volts(1)*(Resist+.65)-V50(1)*.65)/R2 
6560 IF N=38 THEN Pow(I)=(V50(I)-Volts(1))*(Volts(I)*(Resist+.3)-V50(1)*.3)/R2 
6570 ! 

6580 Resist=1.852 

6590 R2=Resist*Resist 

6600 ! 

6610 IF N=24 THEN Pow(1)=(V49(I)-Volts(I))*(Volts(I)*(Resist+. 1)-V49(1I)*. 1)/R2 
6620 IF N=33 THEN Pow(I)=(V50(1)-Volts(I))*(Volts(1)*(Resist+.7)-V50(1)*.7)/R2 
6630 ! 

6640 Resist=1.8333 

6650 R2=Resist*Resist 

6660 ! 

6670 IF N=11 THEN Pow(I)=(V47(I)-Volts(I))*(Volts(1)*(Resist+.5)-V47(1)*.5)/R2 
6680 IF N=12 THEN Pow(I)=(V47(I)-Volts(1))*(Volts(I)*(Resist+. 1)-V47(I)*. 1)/R2 
6690 IF N=13 THEN Pow(1)=(V47(I)-Volts(1))*(Volts(1)*(Resist+.2)-V47(1)*.2)/R2 
6700 IF N=41 THEN Pow(1)=(V51(I)-Volts(I))*(Volts(I)*(Resist+. 1)-V51(1)*.1)/R2 
6710 IF N=45 THEN Pow(1)=(V51(I)-Volts(1))*(Volts(I)*(Resist+.6)-V51(1)*.6)/R2 
6720 ! 

6730 Resist=1.8182 

6740 R2=Resist* Resist 

bio0)! 

6760 IF N=21 THEN Pow(I)=(V48(I)-Volts(1))*(Volts(I)*(Resist+.7)-V48(1)*.7)/R2 
6770 ! 

6780 Resist=1.8 

6790 R2=Resist*Resist 

6800 ! 

6810 IF N=8 THEN Pow(I)=(V46(I)-Volts(I))*(Volts(I)*(Resist+.2)-V46(1)*.2)/R2 
6820 IF N=19 THEN Pow(I)=(V48(1)-Volts(I))*(Volts(1)*(Resist+.2)-V48(1)*.2)/R2 
6830 IF N=20 THEN Pow(1)=(V48(I)-Volts(1))*(Volts(1)*(Resist+.7)-V48(1)*.7)/R2 
6840 IF N=22 THEN Pow(I)=(V48(I)-Volts(1))*(Volts(1)*(Resist+.7)-V48(I)*.7)/R2 
6850 ! 

6860 Resist=1.6667 

6870 R2=Resist* Resist 

6880 ! 

6890 IF N=1 THEN Pow(I)=(V46(1)-Volts(I))*(Volts(1)*(Resist+.3)-V46(1)*.3)/R2 
6900 IF N=3 THEN Pow(I)=(V46(I)-Volts(1))*(Volts(1)*(Resist+.7)-V46(1)*.7)/R2 
6910 IF N=7 THEN Pow(]=(V46(I)-Volts(1))*(Volts(1)*(Resist+.2)-V46(1)*.2)/R2 
6920 IF N=16 THEN Pow(I)=(V48(I)-Volts(1))*(Volts(1)*(Resist+.4)-V48(1)*.4)/R2 
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6930 
6940 
6950 
6960 
6970 
6980 
6990 
7000 
7010 
7020 
7030 
7040 
7050 
7060 
7070 
7080 
7090 
7100 
7110 
ZO 
7130 
7140 
7150 
7160 
7170 
7180 
fi Lek, 
7200 
e210 
(220 
7230 
7240 
7250 
7260 
T2700 
7280 
Zoo 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 
7410 
7420 


IF N=17 THEN Pow(I)=(V48(1)-Volts(1))*(Volts(I)*(Resist+.4)-V48(1)*.4)/R2 
IF N=26 THEN Pow(])=(V49(I)-Volts(1))*(V olts(1)*(Resist+.2)-V49(1)*.2)/R2 
! 


! 
NEXT I 
PRINTER IS 9 
A$="TIME (MIN)" 
B$="TEMPERATURE C” 
C$="POWER WATT" 
Ds= TIVE (SEC) 
E$="NUMBER' 
IF N=) THEN 
PRINT "THIS THERMOCOUPLE IS BAD " 
RINT 
END IF 


IF N=3 THEN 
PRINT "THIS THERMOCOUPLE IS BAD " 
PRINT "" 

END IF 


IF N=6 THEN 
PRINT "THIS THERMOCOUPLE IS BAD " 
PRINT "” 

END IF 

! 

IF N=22 THEN 
PRINT "THIS THERMOCOUPLE IS BAD " 
PRINT ”” 

END IF 

! 


i N=29 THEN 
PRINT "THIS THERMOCOUPLE IS BAD " 
PRINT 

END IF 


PE N=37 THEN 
PRINT "THIS THERMOCOUPLE IS BAD " 
PRINT ** 
END IF 
PRINT "TIME ",TIMES(TIMEDATE) 
PRINT "DATE ",DATES(TIMEDATE) 
ERIN: 
PRINT "HEATER NUMBER °;N 
PRINT ° 
PRINT "CENTER COLUMN PARAMETERS" 
PRINT "CYCLE TIME ",Np 
PRINT "MAXIMUM VOLTAGE |",Maxv 
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7430 
7440 
7450 
7460 
7470 
7480 
7490 
7500 
7510 
7520 
7530 
7540 
7590 
7560 


PRINT "MINIMUM VOLTAGE" Minv 
PRINT" " 

PRINT "SIDE COLUMN PARATERMETERS" 
PRINT "OUTPUT" Scp,'"WATTS" 


PRINT. > 

PRINT "Tamb =",Tamb 

PRINT" " 

PRINT USING "6A,2X,10A,4X,10A,4X,14A,2X,10A";E$,D$,A$,B$,C$ 


Runtime=Tf-Ti 

FOR I=1 TOM 

Tsec(1)=(//M)* Runtime 

Tyme(1)=(1*Runtime/M )/60 

PRINT USING "1X,4Z,2X,7D.D,5X,4D.DD,11X,DD.DD,10X,DD.2D";1,Tsec(I), Tyme(I) 


,T(1)-Tamb,Pow(I) 


7570 
7580 
7590 
7600 
7610 
7620 
7630 


WAIT 0. 

NEXT I 

BEEP 3011.06,1. 
BEEP 406.90,.3 
BEEP 3011.06,.3 
BEEP 406.90,.3 
END 
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202 
254 
260 
262 
290 
300 
310 
330 
340 
350 


APPENDIX D: PROGRAMS FOR THE HP 85 COMPUTER 


TOUTE PPPOE PEPE Peer ee eer e reer etree 


eenereceececeorvreeseeeteaeeeeseeseeeee eee eee nees 


SQURE WAVE 
! TRANSIENT POWER PROGRAM 


THPOPPT PPPOE Otero ieee ee eres 


eeeeeerere ee eee eoeeeeeseeseeeeeevoeeeeeesee 


DIM V1(1399) 

DISP "INPUT MAX VOLTAGE" 
DISP "AND MIN VOLTAGE" 
INPUT R1, R2 

DISP "INPUT CYCLE TIME" 
PNPUT P 

D=.05 

PeP* 1.13 

Z=P/D 

N=INT(Z) 

K=INT(N/2) 

FOR I=1 TO K 
V1I(D=INT(R1*100)+2000 
NEXT I 

FOR I=K TO N 
V1I(D=INT(R2*100)+2000 
NEXT I 

BEEP 

FOR I=1 TON 

OUTPUT 706 USING "#,4D" ; V1(D) 
NEXT I 

GOTO 290 

END 


243 


TEPPER PTEP OPP TT ETE PeR Peer eee eeey 


eeoer eee eee set ene aneer eee eee eevee eee ee eaee 


TRIANGULAR WAVE 
! TRANSIENT POWER PROGRAM 


t 
TRTPPP PEPE EPP R REET EAPO rP rere erey 


eevee e eee ee see eee eee weeesnerseeeereeeeerenes 


DIM (E(1399) 

DIM V1(1399) 

DISP "INPUT MAX VOLTAGE" 
DISP "AND MIN VOLTAGE" 
INPUT R1, R2 

DISP "INPUT CYCLE TIME" 
INPUT P 

D=.05 

| et) Pas Ea Bs 

Z=";D 

N=INT(Z) 

K=INT(N/2) 
V1(K)=INT(R1*10)*10+2000 
V1(0)=INT(R2*10)*10+2000 
E(K)=R1 

E(0)=R2 

L=K-1 

FOR I=1 TO K 
E(1)=(R1-R2)/K+E(I-1) 
E(N-1)=E(I) 
V1I(D=INT(E()*100)+2000 
VI(N-D=V1(1) 

NEXT I 

M=N-1 

BEEP 

FOR I=0 TO M 

OUTPUT 706 USING "#,4D" ; V1(1) 
NEXT I 

GOTO 290 

END 
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350 
340 
350 


CONSTANT 
POWER SUPPLY PROGRAM 


DISP "ENTER THE POWER SUPPLY ’S" 

DISP "FULL SCALE OUTPUT VOLTAGE 

INPUT F 

IF F<0 THEN BEEF ISP "NEGATIVE VOLTAGE’ @ GOTO 60 


DISP "ENTER THE DESIRED OUTPUT VOLTAGE" 

INPUT W 

IF V<0 THEN BEEP ISP "NEGATIVE VOLTAGE NOT ALLOWED" @ GOTO 110 
] 


! CALUCLATE PERCENT OF 

! POWER SUPPLY FULL SCALE 

P=V/F*100 

IF P>99.9 THEN BEEP ISP "VOLTAGE TOO LARGE" OTO 110 
! 

! SELECT RANGE 

IF P>9.99 THEN GOTO 280 

! 


' LOW RANGE CALCULATIONS 
M=INT(P*100+.5) 

D=M+1000 

GOTO 320 

! 


' HIGH RANGE CALCULATIONS 
MalLN F(P*10+.5) 

D=M+2000 

t 


OUTPUT 706 USING "#,4D" ; D 

DISP "THE PROGRAMMED DATA WORD IS";D 
GOTO 110 

END 

GOTO 290 

END 
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APPENDIX E: SAMPLE UNCERTAINTY CALCULATIONS 


Sample calculations are for the top heater location triangular wave power 
pattern baseline case; that is mean power input of 2.0 watts, amplitude 1.0 
watts, and a frequency of pulsation of 0.050 Hz. The uncertainty is calculated 
by considering timewise jitter and bias components for each appropriate 


variable. 


No = Q conv’ 

k A(AT) 

aes 5 @ ane 
‘nee 








Nu _ || 2Qeonw | , (BLY , (BAY , (AT) 
Nu Q) onv L A aE 


1 


= og i BQ) .onv 
6Qeonv = | 50 «| be 204 


(3QY + (SQcona)?l 








| 


OAT) rw = “ sr) sg 220. 
TW 





oT oT 


amb 


(ST + OT, 
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where: 


Nu 


Q 
Glare 
Qeona 
SQpias 
SQrw, 
SQeond 
L 

éL 

SA 

A 

AT 
OAT ras 
omy 


(OT ae) Twa 


1.888 


2.0 Watts 
1.93 Watts 
0.07 Watts 
0.05 Watts 
0.05 Watts 
0.0035 Watts 
2.940 x 10° m 
0.0 

0.0 

186.4 x 10° m? 
20.04 K 

Ok 

On025- 16 


0.025 K 
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8Qoom» = [(0.05)2_45 + (0.05)%pyy + (0.0035)2)" 


0.0708 watts 


(AT) ry, = (0.025)? + (0.025)? 


0.0354 K 
SAT pag = 0.1 K 


éNu _ |(.0708F , 99. 90 « (-0384) . (_-1 
Nu 1.93 20.4 50.4 


SNu 




















= 0.0371 
Nu 
SNu = 0.070 
Nu = 1.888 + 0.070 
Ch gBQL* 
RAV? 
8Gr _|(SQY ,|8L*| , (84 
Gr Q lz A 
where: 
er = 6808 
Q = 2.0 watts 
OGG = 0.05 watts 
OQ my; = 0.05 watts 
L? = 74.71 x 10}? m‘ 
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Si = 0.0 








i SG 10 me 
sy (08) 00a 
Gr 7 2) pws 
oGr _ 9.0354 
Gr 
SGr = 240.7 
Gr = 6808 + 240.7 
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